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BALLISTIC  RESEABCH  LABORATORIES 
THE  PENNSYLVANIA  STATE  UNIVERSITY 
PROJECT  NO:  DA-56-054-0RD-5576-RD 

Progressing  Wave  Analysis  of  Blast  Waves  in  Spheres 
Report  No.  1 

ABSTRACT 

A  theory  of  crater  fomation  hy  inpact  awaits  better  xmderstandlng 
of  the  process  of  shock^wave  propagation  in  solid-: «  especially  waves  wl-th 
spherical  synnetry.  This  paper  studies  theoretically  and  experimentally 
the  non-steady  motion  of  metallic  spheres  initiated  by  e:qploslve  blast  in 
a  spherical  cavity.  Ibe  method  of  progressing  waves  is  applied  to  deter¬ 
mine  the  radius  versus  time  diagram  of  the  propegation  of  the  wave  into 
the  material  and  leads  to  values  for  cavity  sizes.  The  assumptions  are 
made  that  the  material  in  the  vicinity  of  the  ca-^rlty  possesses  a  poly¬ 
tropic  equation  of  state.,  that  entrc^y  is  constauit  for  an  element  of  mater¬ 
ial,  and  that  the  toteuL  energy  is  constant  in  xime.  The  original  inurtlail 
dlfferentleCL  equations  of  the  problem  are  then  reducible  to  a  succession 
of  ordinary  differential  equations.  Using  the  Ranklne-Hugoniot  relations 
as  initial  condltlcms  at  the  shock  front,  these  equations  have  been  inte¬ 
grated  using  a  numerical  progrsoi  developed  for  a  number  of  metals  and  the 
construction  of  r,t -diagrams  ceunried  out  with  -values  for  pairtlcle  -veloci¬ 
ties  and  pressure  -variation  on  -the  inner  surface.  Ihe  solutions  of  -the 
differentied  equation  of  progressing  waves  have  been  studied  by  construct¬ 
ing  a  solution  diagram,  which  is  similar  to  a  hodograph  plane.  An  euaalysls 
was  made  of  this  plane  euid  the  role  played  by  the  singularities  of  the 
differential  equation.  The  appropriate  solution  curve  starts  very  close 
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to  a  positive  node  of  the  equation,  then  approaches  very  close  to  a  saddle - 
type  singularity  in  a  corner  of  the  plane  which  has  been  fovind  to  represent 
a  meaningful  physical  boundary  condition,  namely  that  pressure  emd  velocity 
tend  asyn^jtoti.cally  to  zero  vith  increasing  time  in  the  proper  manner. 

Using  an  equation  of  state  for  edumintm  (obtained  from  data  from  Los 
Alamos  publications)  there  is  obtained  the  expansion  of  a  cavity  of  1.?  can 
radius  in  a  thick  alminum  sphere,  filled  with  31„6  grams  of  Pentolite,  to 
its  final  measured  value  of  J.O  cm.  The  initial  pressure,  which  is  of  the 
order  of  300  kilobars,  drops  to  less  than  100  kb  in  less  than  2  microsec¬ 
onds.  At  this  time  the  shock  velocity  drops  to  its  acoustic  or  elastic 
value  in  the  material.  Hcjwever,  the  cavity  continues  to  expand  to  its  final 
stage  in  a  time  of  80  microsec ..  The  simultaneous  drop  of  pressure,  velocity, 
and  departxire  of  the  medium  from  the  polytropic  equation  of  state  signals  the 
termination  of  the  shock  regime  This  appears  as  a  triangular  region  on  the 
r,t-diagram,  boionded  by  the  shock  front,  the  inner  cavity,  and  the  line  t  =  2 
microsec  Beyond  this  time  the  material  of  the  zone  continues  to  flow  radially 
outwards  essentially  as  an  incompressible  fluid- 

The  assunptlon  that  metals  behave  similar  to  gases  as  a  result  of  an  ex¬ 
plosion  or  ijipact,  is  limited  to  this  shock  zone,  which  is  shorter  in  duration 
than  may  have  been  expected  Prom  the  known  outward  displacements  of  the  out¬ 
er  radius  of  the  sphere  »  1  Inch)  the  increase  in  cavity  volume  is  account¬ 

ed  for  gecmietricaliy 

A  similar  time  scale  of  events  may  be  expected  to  tedie  place  in  eui  im¬ 
pact  crater,  that  is,  the  shock  wave  regime  should  be  terminated  essentially 
before  the  flow  of  material  both  radial  and  tangential,  has  started. 

The  Ballistic  Research  Laboratory,  Aberdeen  Proving  Grounds,  is  conduct¬ 
ing  a  parallel  experimental  program  on  thick  aluminum  spheres  and  have  fur¬ 
nished  data  needed  for  this  analysis 
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1..  INTRODUCTION  -  Description  of  the  Problem 

The  Ballistic  Research  Laboratory  (BRL),  Aberdeen  Proving  Grounds, 

Is  engaged  In  conducting  an  experimental  program  of  Internal  explosions 
in  small  cavities  In  metal  spheres.  The  aim  of  this  report  Is  to  give 
the  restilts  of  an  analytical  study  of  the  problem,  with  particular  em¬ 
phasis  on  the  propagation  of  shock  waves  in  the  metal  Immediately  after 
detonation  of  the  explosive. 

Figs.  1  and  2  show  a  7  inch  sphere  of  aliunlnum  such  as  has  been  used 
to  date  In  BRL  experiments.  The  phenomena,  and  hence  their  analysis, 
which  result  from  detonation  of  the  explosive  In  the  inner  cavity,  are 
cotnplex  because  different  effects  predominate  In  different  parts  of  the 
material  Thus,  there  is  an  Innermost  zone  or  spherical  shell  where  very 
large  radial  displacements  have  occurred  under  temperatures  and  pressures 
far  beyond  the  range  of  conventional  mechanical  behavior.  The  material 
Is  In  some  type  of  "fluid"  state  in  this  zone  and  shows  relatively  little 
tendency  for  cracks  to  Initiate  there.  Next,  there  Is  em  Intermediate 
zone  of  the  sphere  where  dlspXacements  have  dropped  to  elastic  ranges  and 
where  brittleness  seems  to  have  returned,  as  evidenced  by  the  many  small 
tension  and  shear  cracks  which  have  formed  there.  A  few  of  the  stronger 
cracks  which  get  started  may  penetrate  Into  the  adjacent  regions  and  ulti¬ 
mately  reach  the  boundaries  Finally,  there  Is  an  outermost  zone  dominat¬ 
ed  by  the  effect  of  the  external  boundary  of  the  sphere.  Here  reflection 
effects  such  as  scabbing  cracks  are  often  observed. 

The  transition  between  successive  zones  are  not  exact,  but  in  some 
specimens,  rather  surprisingly  enough,  are  fairly  sharply  delineated. 

In  this  report  we  shall  make  an  analysis  of  the  Innermost  shock  zone. 
Its  direct  aim  is  to  provide  a  description  of  the  shock  process  In  the 
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metal.  More  specifically  a  useful  theory  must  furnish  a  time  for  the 
duration  of  the  process,  the  size  of  the  zone  influenced  hy  the  shock 
front,  values  for  the  displacements,  and  thezmoflynaffllc  variables  of  pres¬ 
sure,  density,  and  tenqperature  in  the  material.  A  useful  tool  is  the  r,t- 
dlagram  vhlch  shove  the  path  made  by  a  set  of  concentric  spherical  shells. 
This  diagram  is  possible  because  of  the  single  space  coordinate. 

The  problem  of  the  shock  e3q>ansion  of  spherical  cavities  is  closely 
related  to  that  of  crater  fozmtion  by  hypervelocity  projectiles.  The 
features  ve  have  outlined  above  are  present  in  the  crater  problem  as  veil. 
The  crater  problem  carries  vlth  it,  however,  the  further  cospllcation  of 
tangentlail  flow,  thus  requiring  two  space  coordinates.  Sxc^t  for  the 
presence  of  the  plug  shown  in  Figs.  1  and  2,  the  arramgement  for  the 
blasts  have  spherical  symmetry,  and  ve  amy  confidently  assert  that  radial 
motion  occurs,  so  that  aiU  the  physicail  quantities  of  the  problem  depend 
on  only  one  space  coordinate.  There  is  however,  a  slight  actual  departure 
because  of  the  plug  or  because  of  asymmetrical  detonation,  and  which  are 
not  Important  to  the  problem. 

2.  Spherical  Shocks  -  Time  Sequence  of  Effects 

Just  as  the  study  of  the  problem  has  been  conveniently  divided  up  in¬ 
to  spatlsd  zones,  ve  can  divide  up  the  sequence  of  events  in  the  spherical 
blast  process  for  detailed  sinsilysls  as  follows: 

a)  Initial  Stage  -  Here  the  detonation  wave  of  the  ejqplodlng  gas 
starts  contact  vlth  the  solid  and  then  generates  a  shock  wave  in  the  solid. 
This  stage  mic^t  be  considered  as  terminated  when  the  density  in  the  solid 
has  dropped  to  its  free  space  value,  at  the  inner  cavity. 

b)  Eicpansion  Stage  -  The  coDq)reBsed  solid  e^qiands  radially  outward 

and  actually  forms  the  cavity.  This  stage  is  dominated  mostly  by  inertia 
forces . 


-5- 


c)  Final  Stage  *  Here  the  shocit  va.ve  decays.,  permanent  deformation 
of  the  cavity  stops,  and  the  material  has  undergone  some  permanent  plastic 
strain. 

The  first  stage  lasts  up  to  about  ,2.5  micro-seconds.  The  second 
turns  out  to  he  relatively  long  and  can  take  up  to  about  100  microseconds 
or  even  longer.  It  must^  of  course,  be  understood  that  these  phases  need 
not  be  distinctly  separated  events  in  time,  especially  the  terminating 
phase  of  the  expansion. 

5.  Basic  Mathematical  Equations  of  Shock  Waves 

Because  of  spherical  symmetry  our  problem  is  reducible  to  a  radial 
and  a  time  coordinate.  Shock  wave  propagation  In  a  solid  Is  very  closely 
related  to  that  of  a  spherical  wave  In  a  gas.  We  nay  make  the  following 
assunqptlons  about  the  medium; 

1.  Thermodynamic  eq,uillbrlum  holds  (See  [2],  p.  3)>  l»e.  that  changes 
of  state  are  adiabatic  By  this  we  mean  that  entropy  Is  constant  along  a 
"particle  path",  i.e  ,  a  fixed  element  of  the  medium. 

2.  The  medium  i£  a  perfect  fluid,  1  e  ,  any  rigidity  or  shear  effects 
are  neglected. 

5-  The  effects  of  entropy  changes  are  negligible,  i.e.,  that  the 
pressure  is  a  ranotion  of  the  density  alone  (the  medium  is  said  to  be 
barotroplc ) . 

The  total  energy  available  for  the  motion  is  fixed. 

The  conservation  Laws  for  an  element  of  material  expressed  in  Lagrangean 
form,  i.e.,  along  the  p)artlcie  paths,  are  as  follows: 
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dp/dt  =  -  p(du/5r)  -  2up/r  s 

-  (p/r^)d(r^u)/^ 

(mass) 

(3.1) 

du/dt  =  -  bp/phr 

(momentTim) 

(3.2) 

f(P,p)  »  0  or  df/dt  -  0 

(eq.  of  state) 

(3.3) 

Thus  If  the  medium  were  assumed  to  he  polytrpplc  with  the  adiabatic 

ent  7,  we  would  have,  for  (3*3)> 

esqpon- 

f(p»p)  ■  PP  ^  ■  const.  «  A 

(3.4) 

These  equations.  In  Eulerlan  form, 

tlves,  become 

with  subscripts  denoting  partial  derlva' 

p^  +  u  p^  +  pu^  +  2up/r  -  0 

(mass) 

(3.5) 

u^  +  uUy  +  p^p  •  0 

(momentum) 

(3.6) 

(pp'^')^  +  u(pp"’')^  -  0 

(state) 

(3.7) 

The  third  of  these  eqiiatlons  Is  not  quite  equivalent  to  since  It 

only  expresses  the  fact  that  the  entropy  Is  constant  along  the  path  of  an 
element,  euid  does  not  lnq>ly  Its  constemcy  throu^out.  This  Is  a  difference 
from  the  case  of  plane  waves;  euiother  difference  from  the  equations  of  one- 
dlmenslonal  flow  is  the  additional  term  2up/r  occurring  In  (3*1)  and  (3- 5)/ 
which  stands  essentially  for  the  spherical  attenuation  of  the  wave.  This 
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term  of  cotirse  Is'  very  important  to  the  problem. 

A  conplete  geometrical  description  of  the  disturbance  is  afforded 
by  the  construction  of  an  r,t-diagram,  as  shown  in  Fig.  ?•  Here  the  solid 
lines  represent  the  motion  of  the  points  of  a  spherical  surface,  referred 
to  as  a  "particle".  The  most  prominent  feature  in  this  diagram  is  a  dis¬ 
continuity,  or  shock  front  which  propagates  through  the  materieJ.  at  the 
head  of  the  disturbance.  This  curve,  together  with  the  cavity  boxuidary, 
defines  a  region  (shaded  in  the  figure)  in  which  the  solution  to  the  system 
of  partial  differential  equations  (5*5)>  (3.7)  applies.  Certain  boundary 
conditions,  to  be  discussed  later,  must  be  satisfied.  However,  the  diffi¬ 
culty  of  the  problem  is  that,  unlike  the  conventional  boundary  value  prob¬ 
lems,  here  the  boundary  curves  are  themselves  unknown,  and  must  be  found 
as  part  of  the  problem.  In  fact,  the  determination  of  these  two  curves 
are  the  most  In^jortant  part  of  the  problem. 

k.  The  Method  of  Progressing  Waves 

The  idea  of  this  and  similar  mathematlcsil  methods  is  to  reduce  the 
partial  differential  equations  to  ordinary  ones.  This  is  accooiplished  by 
assuming  the  specific  form  for  the  shock  front  curve  and  imbedding  it  in 
a  one -parameter  family  of  curves.  These  curves  are  called  "progressing 
waves".  For  general  details  of  the  method,  see  [2]  p.  kl9-J«-33.  The 
method  was  used  by  R.  G.  Newton  [h]  to  analyze  blast  shock  problems. 

Our  "progressing  wave"  solutions  are  defined  to  be  of  the  form, 

u  =  t^  I  U  (|)  (4.1a) 

with  I  =  rt”® 

P  =  D  (6)  (4.1b) 

p/p  =  t®  5^  P  (0 


(4.1c) 


-8- 


where  a,  p,  8,  e  are  parameters,  eund  U,  D,  P  functions  to  "be  deter¬ 
mined.  By  introducing  this  variable  t  ve  have  defined  geometrically  a 
family  of  surfaces  5  =  const,  in  the  r,t -plane,  which  will  play  an  im¬ 
portant  role  in  the  analysis.  Although  these  are  not  the  trajectories  of 
the  particles  of  the  medium,  we  shall  see  that  the  shock  front  belongs  to 
this  family  of  siirfaces. 

Ve  now  explore  these  solutions  mathematically  by  substituting  the 
expressions  (k.l)  into  the  equations  of  motion  (3.5)-(3-7)>  giving  respect¬ 
ively. 


5t^‘^  j^pU  -  a  (5U'  +  U)  +  t^*^^  U(tU'  +  U)  + 

(2P  +  J  P  D'/D  +  i  P')  ]  =  0  (4.2a) 

t®‘  ^  t(U  -  a)  i  D'  +  8  D  +  (j  U’  +  3U)D]  »  0  (4.2b) 

^8+€  -  1  ^^2  p  jj  _  ^  jjp 

-  a(2DP  +  J  D'P  +  {  DP')]  +  t^"°‘'*'^[-  PD' 

+  (  (2DP  +  5  D'P  +  6  DP')]  t  U  ]■  .  0  (4.2c) 


The  sense  in  these  equations  is  that  it  is  possible  to  eliminate  the  explicit 
fiactor  t  by  properly  choosing  the  esqponents,  thereby  leaving  a  system  of 
functions  of  one  independent  variable  (.  This  is  accomplished  by  letting 


c-2p  ;  P-a-1 


(4.3) 
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80  that  ^  ^^^  =  1  and  we  then  have  (after  also 

dividing  by  Jt^  t^  t^  ^  ^  ^  which  are  not  zero  for  t  >  0, 

and  seme  simplification, ) 

PU  -  a  (|  U'  +  U)  +  U(S  U'  +  U)  +  (2P  +  5  PD'/D  +  e  P')  =  0  (1^.4a) 

(U  -  a)  I  D'  +  6  D  +  (|  U*  +  3U)D  =  0  (4.4b) 

P'  j(U  -  a)  +  p[2p  -  8(7  -  1)  +  2(U  -  a)]-(D'/D)p  s  (7  -  1) 

(U  -  a)  =  0  (4.4c) 

We  now  have  a  system  of  ordinary  differential  equations  for  the  \mknown 
functions  U({),  0(5),  P(|)  suid  two  free  parameters  a  and  8.  The  sub¬ 
stitutions  (4.1),  which  may  appear  artificial,  are  thus  Justified. 

Ve  shall  now  reduce  the  nTsnber  of  variables  further.  Solving  (4.4b) 
for  5  D'/b  gives 

(  D'/D  =  -  (8  +  I  U'  +  3U)/(U  -  a)  (4.5) 

When  this  is  put  into  the  remaining  two  equations  we  obtain  the  pair, 

(u  -  a)  (t  u'  +  u)  +  p  u  +  P  [2  -  (8  +  5  u'  +  5U)/(u  -  a)] 

+  I  P’  =  0  (4.6a) 

-  (r  -  1),  (8  +  5U”+3U)-2p+(7-l)6.2(U-a)- 

-  (U  -  a)  j  P'/i*  =  0  (4.6b) 
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These  equations,  linear  in  JU'  and  JP'  may  be  simultaneously  solved, 
giving 


I  U'  =  [  -  U  (U  -a)  (U  -  1)  +  (B  +  2P  +  3U  7)P]A  (^-Ta) 

S  P'  =  P  |(7  -  1)U(U  -  1)  -  (U  -0!)[(37  -  1)U  -  2]  + 

+  ~  +  2r)P  ]  /A  (J^.Tb) 

where  A  b  (U  -  O.)^  -  7P  (^'Tc) 

The  final  step  is  to  obtain  by  division  of  (4.7b)  by  (4.7a)  the  ordinary 
differential  equation 


dU 


(4.8a) 


where,  after  simplification, 

»(U)  »  7  U(3a  -  1  -  2U)  +  (3  -  ot)  U  -  2a 

Q  (u)  -  {[2P  -  (7  -  1)  Bl/(u  -  a)}  +27 

^  (4.8b) 

R  (U)  -  U(U  -  a)  (1  -  U) 

S  (U)  -  6  +  2P  +  3U7 


This  is  the  basic  differential  equation  for  progressing  waves.  After  the 
appropriate  solution  has  been  foxmd  for  P  =  P(tJ),  the  function  |  »  |(u) 
is  found  by  a  quadrature  of  (4  7a)  and  the  density  function  D({),  from 

(4.5). 
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These  progressing  vave  solutions,  as  we  shall  see,  provide  a  sufficient¬ 
ly  general  mathematical  description  of  an  expanding  cavity  reasonably  consis¬ 
tent  with  the  given  conditions  of  initiation  of  the  process.  There  ronalns 
the  problem  of  choosing  the  two  parameters  0!  and  6.  For  this  we  have 
two  possibilities  of  an  assxm^tlon: 

(a)  the  motion  is  Isentroplc 

(b)  the  motion  is  adiabatic,  with  constant  total  energy. 

If  (a)  holds,  then  from  eqs.  (3«^)  and  (4.1), 

pp"’’  -  (t*{^  P)  A 

or  t'  *  -’■)»  A 

requiring,  for  Independence  of  time,  that 

€  +  (1  -7)6  -  0.  (4.9) 


This  condition  is  not,  in  general,  satisfied  by  a  spherical  wave.  Instead, 
we  have  the  condition  (b),  that  is,  constant  total  energy,  ^is  is  a  reason¬ 
able  one  for  the  cavity  expansion  process,  because  of  its  short  duration, 
provided  certain  secondary  effects  are  neglected.  With  (  ■  represent¬ 
ing  the  shock  front  at  a  time  t,  the  total  energy  in  the  fluid  shell  (poten¬ 
tial  +  kinetic)  at  time  t  is  given  by 


(4.10) 


a  ,  , 

where  r^  =x  I  t  is  the  inner  radius  of  the  shell  (Fig.  1)  and  r,  ■ 

a 

(^t  is  the  location  of  the  shock  front. 
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The  first  term  arises  from  the  polytropic  relation  p  ■  Ap^  which 
has  been  assumed  for  the  material  and  the  fact  that  the  work  done  by  the 
material  in  becoming  compressed  from  an  initial  volume  to  a  volume 
V  is  given  by 

V  P 

E(V)  «  J  pdV  -  p^V^  J  A  p^(  -  const. 

"o  **  '  (k.n) 

where  we  have  used  the  relation  pV  Using  the  substitutions  In 

(k.l)  the  energy  e^gpression  becomes 

E(t)  ~  f  ^(t®  ®  ^  +  I  t®  U^)5^(6)r^dr 

r 

o 

Substituting  r  ■  {t®*,  dr  ■  t^dj,  since  t  is  constant, 

E(t)  -  k«t®  5a  -  2  J  +  i  u2)  D(5)  (4.12) 

*0  ■ 

Since  the  Integral  is  Independent  of  t,  we  make  the  energy  independent 
of  time  by  satisfying  the  relation 


6  +  5a  -  2  -  0 


or 


6  -  2  -  5a 


(^^•13) 
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5.  Boundary  Conditions  at  Shock  Front 

We  shall  narrow  down  the  nvuriber  of  parameters  hy  examining  the  com¬ 
patibility  of  our  solution  with  the  basic  Rankine -Hugoniot  conditions 
across  a  shock  front.  If  the  undlsttirbed  and  disturbed  medium  parameters 
are  u^,  p^,  p^  and  u^,  p^  respectively  and  the  shock  wave 

velocity  is  C,  then  these  relations  are  (See  [2],  p.  123-4), 

p^  (u^  -  C)  =  p^  (uj^  -  C)  =  m  (conservation  of  mass)  (5 -la) 

p,  u,  (u,  -  C)  -  p  u  (u  -  C)  =  p  -  p,  (conservation  (5*lb) 

^  ^  ^  °  ^  of  «c«ent>M) 

“l  <  2  ("l  •  ■  ‘’o  <  f  <"o  -  =>  •  Po  “o  • 

-  p^  u^  (conservation  of  energy)  (5-lc) 

where  E  =  for  a  polytrpplc  medium. 

7-1  p 

When  the  undisturbed  state  is  a  medium  at  rest,  with  u^  s  0,  these 
reduce  to 

Pi  (C  -  Uj^)  -  p^C  «  0  (5-aa) 

Pi  (C  -  u^)  -  (Pj^  -  Pq)  =  0  (5.2b) 

Pi  (  I  ^1^  +  ^  )  (C  -  u^)  -  p^  u^  =  0  (E^  =  0)  (5.2c) 

It  is  useful  to  have  (5.2a)  and  (5.2b)  solved  for  the  velocities,  giving 


=  (Pi  -  V^)Ap^C) 


(5.5a) 
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(5.3b) 


If  p  is  in  kilobars,  p  in  gm/cc,  then  u  and  C  will  ccame  out  in 
km/sec  by  the  following  formulas: 


Ui  -  -  Po)/(10  pQ  C) 


(5.i4a) 


(5.4b) 


It  is  useful  to  have  these  relations  for  a  polytrppic  medim 


(5.5a) 

(5.5b) 


From  (5.3a),  with  p  »  y  (7-l}/(7+l)  ,  if  p^  is  negligible. 


u^  -  C  (1  -  V^) 

(5.6a) 

Pi  -  p^  u^  C  -  p^  (1  -  /) 

(5.6b) 

Pl/Po  =  (7  +  l)/(7  -  1)  =  l/p^ 

(5.6c) 

c  =  ^  dp/dp  =  C 

The  last  quantity  is  conveniently  referred  to  as  the  "so\ind  speed"  in 
shock  wave  emalysis. 
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These  relations  apply  Just  as  well  to  a  spherical  or  curved  surface  as  to 
a  plBuie,  since  the  effect  of  spherical  divergence  (the  2u/r  term)  on  a 
finite  or  sudden  Jump  is  of  higher  order.  This  may  also  be  shown  geomet¬ 
rically  by  considering  an  infinitesimal  surface  element  of  the  shock  front. 
Since  =  rt  along  the  shock  front, 

0  =  t"^  dr  -  Ctrt"^  ^dt 

so  c  =  ^  =  0(rt  ^  =  a|^t^”^  .  From  (4.1)  and  the  relations  (4.5)  with 
c  -  u^  =  Ctl^t^  ^-t^  (a  -  U),  we  obtain 

|^D(a  -  U)  -  =0  (5.7a) 

(a  -  U)  -  I^^DP  »  0  (5.7t) 

(  I  ~  P)  (a  -  u)  1^^  DPU  =  0  (5.7c) 

We  note  that  the  time  factor  cancels  in  (5'Tb)  and  (5'7c),  so  that  they 
are ' automatically  satisfied,  but  to  secure  independence  of  time  in  (5.7a), 
it  is  necessary  to  make 

6-0  (5.7d) 

With  this  condition,  and  the  relations  (4.5),  the  assumed  fom  for  the 
progressing  wave  solutions  reduce  to 


-i6- 

|o(s) 

p .  ( 

l)  I)(S)P(t) 

I>(s) 

s  -  ( 

if 

vith  I  ■  rt"^ 


(5.8) 


^nils  solution  shovs  that  on  the  shock  front  or  free  surface^  where  |  Is 
constant,  the  physical  quantities  such  as  velocity,  pressure,  density,  and 
wave  velocity  are  constant  on  the  rays  r/t  >  constant.  This  also  dlmen> 
slonallzes  the  functions  (5 “8)  correctly,  for,  with  |  having  the  dimen¬ 
sions  of  length  per  (time)®,  and 

I  =  [LT"®]  D(|)  =  [ML‘^1  density 

U(|)  -  [1]  P(6)  -  [1] 

then  u,  p,  P,  p/p  come  out  correctly  to  have  respectively  the  dimen¬ 
sion  of  velocity,  pressure  density,  and  velocity  squared.  The 

C0Bq)lete  set  of  e^q^onents  is  now 


a  -  2/5  €  -  -  6/5 

P  -  -  3/5  6-0 


(5.9) 


Initial  Conditions 

Since  t  "  (j.  ^  shock  frcmt,  we  have,  Just  behind  it, 

U(6j^)  -  0(1  -  M^)  -  aa/(7  +  1)  (5.10a) 


D(5i)  -  Pq/p^ 

P(5l)  = 


(5.10b) 


(1  -  n^)  -  20^//i7  +  1) 


(5.10c) 
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with  a  =  2/5  or  l/k  according  as  the  energy  or  momentum  condition  holds. 

The  right  side  of  equations  (5-10)  give  us,  for  a  specified  material, 
a  definite  point  in  the  P-U  plane,  through  which  a  single  solution  curve 
is  determined  in  general.  We  may  refer  to  this  point  as  our  "initial  point" 
and  proceed  to  draw  the  solution  curve.  Note  that  the  constant  still 

undetermined,  is  not  needed  for  this.  We  will  discuss  in  Section  10  how 
this  constant  may  he  determined. 

6.  Equation  of  State  of  Aluminum 

In  Fig.  4  is  shown  the  relationship  between  pressure  and  relative 
density  for  aluminum  in  the  range  between  100  and  VOO  kilobars.  The  three 
numbered  data  points  are  taken  from  Ref.  [3]  as  follows: 

24  ST  eJ.vmlnum 

=  2.785  free  space  density 

Cp  »=  0.23  8p.  heat  at  const,  pressiare 


^  (^1  -  69.0  X  10' Vc 

'  /p 


^2 

Pressure  (kb) : 

355.8 

222.7 

155.5 

Rel.  volume  V/V^: 

0,7874 

0.8555 

0.8696 

Shock  wave  velocity  U  : 
(km/sec)  ® 

7.551 

6.927 

6.500 

Free  surface  velocity: 
(km/sec) 

5.250 

2.519 

1.700 

These  points  are  plotted  on  logarithmic  paper  in  Fig.  4,  and  a  straight 
line  fit  throu^  these  points  is  possible,  leading  to  the  polytropic -type 


relation 
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P 


52.7 


100  ltb<rp<400  Kb 


(6.1) 


vhere  p  Is  in  Kllobars. 

A  few  points  are  shown  beyond  the  400  Kb  range  based  on  additional 
data  taKen  from  [3].  Here  there  Is  a  slight  bxtt  consistent  departure  fraa 
the  equation  of  state  (6.I).  However  under  the  conditions  of  our  explosion 
the  range  of  presstires  do  not  exceed  400  lb. 

Below  about  100  Kb  we  have  a  transition  to  elastic -plastic  or  elastic 
behavior.  The  ziattu«  of  this  transition  Is  considerably  uncertain.  Ve  may 
also  note  that^  unllKe  gases j,  the  value  of  7  Is  very  high,  l.e.,  relative¬ 
ly  small  density  changes  occur  under  very  high  pressures. 

The  polytrqplc-type  equation  of  state  Is  a  very  convenient  one  to  use 
for  metals  provided  the  pressure  range  is  restricted,  such  as  In  (6.I). 

It  has  the  advantage  that  the  progressing  wave  procedure  In  eqs.  (4.2) 
can  be  carried  through.  However,  such  an  eqmtlon  of  state  must  always 
be  modified  at  low  pressures  since  the  density  of  a  solid  does  not  tend 
to  zero  with  the  pressure.  Othel-  equations  of  state  have  been  used,  e.g. 
Sedov  In  [11],  and  StanuKovlch  [131,  bave  used  the  fozvula 


p  -  A  [(p/Pq)’'-  1] 


(6.2) 


For  altsBlnum,  the  values  A  «  I87,  7  •  4.27,  p^  •  2.7  provide  a  good 

fit  to  the  data  points  of  Fig.  4. 

Further  discussion  of  equations  of  state  for  solids  Is  given  by 


Huang  [14]. 
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7.  The  Explosion  Process 

A  physical  description  of  the  explosion  process  taking  place  in  the 
expsmding  inner  cavity  is  natttrally  very  Important  to  our  problem  as  this 
provides  the  pressure  forcing  function  required  to  e^qwuid  the  metal.  Al¬ 
though  much  is  knovm  about  the  unconfined  explosion^  the  main  difficulty 
of  our  problem  is  that  the  gas  is  confined  by  the  metal  daring  the  process. 
Our  knowledge  of  the  effect  of  this  confinement,  at  least  for  the  explos¬ 
ives  used  in  the  experimental  study,  is  limited  to  enqpirical  factors. 

A  condensed  explosive  such  as  "Pentollte”  is  converted,  upon  explosion, 
into  a  gas  at  high  temperature  across  a  rapidly  moving  discontinuous  front, 
referred  to  as  a  detonation  wave.  Behind  this  front  the  chemicsd  ccnposi- 
tlon  changes  until  a  number  of  stable  end-products  are  reached.  The  detona¬ 
tion  process  in  Pentolite  has  been  analyzed  in  detail  by  Sheer  [5],  [6] 
using  the  hydrodynamic  theory  of  detonation,  aud  an  equation  of  state 
worked  out  giving  the  pressure-density  relationship  behind  the  detonation 
front.  This  is  shown  in  Fig.  3^  It  may  be  assxmted  that  the  gas  expands 
Isentroplcally . 

The  detonation  values  for  Pentolite  are,  from  [6]: 


Loading  density  of  solid; 
Initial  e3q>lo8lon  pressure: 
Explosion  tenq>erature : 
Detonation  velocity; 


=  1  655  gm/cc 
P^  =  251.250  kllobars 
-  3367.7°  K 
D  •  7807  m/sec 


The  equation  of  state  c\irve  of  Pig.  5  is  not  fitted  too  well  for  the  en¬ 
tire  reuige  by  a  polytropic  equation  (strai^t  line).  However,  over  a 
small  part  of  the  range  of  pressures,  we  may  refer  to  the  slope  7  as 
the  "local"  adiabatic  exponent.  This  varies  from  7  =  3»29  at  the  front 
down  to  7  =  1.52. 
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Por  rough  approximate  purposes  the  hi^  value  of  7  may  be  used  for 
pressures  down  to  100  kb- 

These  results  are  in  very  good  agreement  with  experiments  for  xmcon- 
fined  explosions,  or  explosions  in  air.  However,  there  is  still  lacking 
a  detailed  study  of  the  interaction  between  the  expanding  gaseous  products 
of  the  e3q)losive  and  the  metal  surface  confining  it.  The  expsuision  of 
the  gas  is  considerably  complicated  by  the  detonation  wave  reflecting  off 
the  metal  boundary,  probably  many  times.  Thus  the  loading  pressure  P 

r 

against  the  metal  will  be  higher  than  the  detonation  pressures  P  behind 

s 

the  front,  as  exhibited  In  Pig.  5-  While  the  latter  are  acc\irately  known, 
the  loeidlng  presstires  are  uncertain.  Dorlng,  in  [10],  p.  226,  formula  (4), 
gives  data  of  this  type  for  other  materials  than  Pentolite.  With  these 
data  as  a  guide  an  estimate  of  26^  Increase  in  pressure  due  to  confine¬ 
ment  is  tentatively  STiggested,  i.e. 

P^  =  1.26  P^ 

r  8 

Thus,  for  a  detonation  pressure  of  2^1  kb, 

P  «  292  kb 
r 

The  scune  reference  also  furnishes  the  value  u  *  1400  m/sec  for  the  initial 
peo^lcle  velocity  in  SLlumlnum. 

In  the  early  stages  of  the  eaqploslon  process  (up  to  approximately  2 
mlcrosec)  the  actual  value  may  be  substantie^Jy  higher. 

It  may  reasonably  be  assumed  that  this  factor  remains  constant  at 
different  pressures,  so  that  the  curve  of  Pig.  5  is  Just  shifted  tg>  a 
constant  amount.  Until  more  detailed  studies  are  available,  this  would  be 
the  easiest  way  to  ex+-rapclate  the  confinement  factor  down  to  lower  press¬ 


ures. 
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The  analysis  of  the  confinement  problem  requires  a  model  for  the  den¬ 
sity  variation  with  radl/us  of  the  exploding  gas.  Thus,  the  simplest  model 
would  be  a  gas  in  equilibrium;,  etnd  hence  at  uniform  density,  In  the  cavity, 
as  shown  by  curve  (a),  Fig  6.  Given  the  pressure -density  relation,  the 
pressure  against  the  metal  would  be  easily  calculated  from  the  Instantaneous 
cavity  radius  R(t)  according  to  the  relation  p  =  {R A  more 
realistic  distribution  of  non-uniform  density  is  shown  in  (b),  which  puts 
the  maximum  density  and  pressure  at  the  interface.  It  assumes  the  detona¬ 
tion  front  remains  in  contact  with  the  interface  and  therefore  the  Hugoniot 
relation  of  Fig.  5  ®ay  be  used.  However,  farther  assumptions  are  required 
to  relate  the  density  to  radius  R(’t.!  or  to  the  time  t  itself.  The 
third  model  (c),  which  is  only  sketched,  is  a  further  refinement  in  that 
it  accounts  for  the  reflection  of  the  detonation  front  off  the  interface. 

The  equilibrium  model  {a)  can  only  bsccme  accurate  during  the  latter 
phase  of  the  cavity  expansion,  after  several  reflections  of  the  detonation 
wave.  Fig.  12  conqaares  gas  pressures  on  this  model  with  that  of  the  solid. 
They  are  too  low  at  first,  and  too  high  later.  The  critical  time  occurs 
at  2.5  microseconds.  The  probably  state  of  affairs  is  that  the  detonation 
wave  remains  in  contact  with  the  mets-l  for  thie  period  of  time  and  then  re- 
leues  the  pressure  . 

8.  The  P,  U-  Diagram 

Using  the  condition  cf  constant  energy,  the  2/5 -power  law  holds  and 
the  differential  equation  yk  6}  for  progressing  waves  may  be  solved.  This 
was  done  numerically  by  means  of  Program  I  given,  in  the  Appendix  A  family 
of  integral  curves  in  the  diagram  are  ehown  in  Fig  7  Ail  of  them  issue 
from  the  singular  point  A  enclosed  in  a  rectangle  on  the  diagram,  and  which 
is  located  by  equating 
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F(U,  P)  =  0  ,  G(U,P)  *  0 

in  (4.8) «  For  our  set  of  constants  this  point  is 

=  0.0916  j,  P^  =  0.02884  (8.1) 

A  more  detailed  ajialysls  of  the  singular  points  is  given  in  the  Appendix. 
The  immediate  neighborhood  of  this  point  is  shewn  plotted  in  Fig.  8.  It 
can  be  seen  that  all  the  curves  come  out  of  this  point  along  a  common  tan¬ 
gent,  and  the  point  Is  labelled  an  "unstable  nodal  point”  or  simply,  a 
"source".  The  point  U  «  0,  P  =  0  is  a  stable  nodal  point,  or  "sink", 
and  the  point 

P^  -  0  ,  U  -  0.4 

o  o 

is  a  saddle  point.  (See  Appendix  C.)  The  solution  curve  for  the  physical 
problem  starts  close  to  the  soxirce,  then  runs  over  very  close  to  the  saddle 
point.  There  is  a  unique  curve  (marked  C  in  Figs.  7  and  8)  vhlch  actually 
runs  into  the  saddle,  as  shown  in  Appendix  C.  The  accuracy  of  the  compu¬ 
tations  is  not  sufficient  to  distinguish  whether  our  solution  actually  co¬ 
incides  with  C  or  not. 

The  "source"  point  Itself  represents  physically  an  infinite  shock 
strength,  as  sketched  in  Fig.  9(b).  Here  the  pressure  and  density  Just 
behind  the  shock  front  are  infinite  (with  a  finite  total  impulse,  however) 
and  the  particle  velocity  is  equal  to  the  shock  wave  velocity.  This  type 
of  condition  arises  in  stress -wave  propagation  problems  as  well,  in  the 
foim  of  a  5 -function  at  the  wave-front.  See  [16].  It  is  a  consequence 
of  the  assumed  Instantaneous  (i.e.  step)  loading  of  the  material.  With 
such  a  loading  we  must  start  the  propagation  of  either  a  zone  of  infinite 

pressure  if  the  velocity  is  finite,  or  our  front  must  start  out  with  an 
infinite  velocity. 


-23- 


The  end-point  of  the  solution  curve  C  at  U  =  a,  P  --  0  provides 
a  very  rea.sonable  physical  condition  of  asymptotic  character.  For,  if  we 
consider  any  point  A  in  the  r,t-plane.  Fig.  11,  che  "particles"  must  cross 
the  family  of  curves 

r  =  It®  (8.2) 

(shown  dotted)  from  left  to  right,  since  we  have  conpression  shock.  Thus 
the  particle  curve  (solid)  has  a  lower  slope  at  A  than  that  of  the  dotted 
curve: 

u  =  £u(s)<||  .c..4f 

This  condition  is  always  satisfied  since  all  the  solution  curves  in  the 
U,P-plane  coming  out  of  the  unstable  nodal  point  lie  to  the  left  of  the 
vertical  line  U  =  0.4.  However,  for  the  curve  C  u  -»  0.4  r/t  as 
t  -»  oO ,  i.e.  the  particle  curve  is  asymptotic  to  <8.2).  The  P  =  P(U) 
curves  which  end  in  U  =  0,  P  =  0  give  particle  curves  which  cross  all 
the  curves  of  (8.2). 

9-  The  Radius -Time  r,t -Diagram 

Figs.  10  and  11  show  an  r,t-diagram  plotted  for  an  initial  cavity 
radius  of  I.698  cm.  Pressures  in  kilobars  are  also  shown  on  Fig.  10. 

This  figure  gives  much  more  detail  of  the  early  phase  of  the  expansion 
up  to  t  =  2.5  raicrosec.  from  its  start  at  t  =  0.9  sec.  In  this  elapsed 
time  of  1.6  ^  sec  the  inner  cavity  has  only  grown  to  1.86  cm,  which  is 
only  15^^  of  its  ultimate  change  However,  the  pressure  has  already  fallen 
considerably.  At  the  cavity  surface  it  is  down  to  50  kb. 

We  also  note  that  the  shock  velocity  at  point  labelled  P  or;  the  dia¬ 
gram  is  equal  to  the  known  elastic  wave  velocity  of  the  material..  Beyond 
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this  point  the  2/5 -power  law  for  the  shock  front  starts  to  deviate  from 
this  velocity.  Such  a  condition  represents  a  discrepancy  of  the  progressing 
wave  method  from  this  point  on  which  is  inevitable  because  of  the  equation 
of  state  used. 

Figs.  12,  13,  and  111  show  how  the  pressure,  particle  velocity  and 
density  decay  with  time  at  the  inner  cavity  surface.  We  note  that  p  = 
occurs  for  t  e  2.5  p  sec.  Of  course,  we  may  not  conclude  that  p  becomes 
less  than  the  free  space  density  beca\ise  the  equation  of  state  (6.I)  no 
longer  applies. 
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10.  Results 


TABLE  1 


Variation  of  physical  quantities  along 
solution  curve  and  cavity  surface 


T 

u 

P 

1 

E 

sec  X  10 

mis 

ergs  X  10 

C/wu 

.8684 

.09302 

.02856 

451.28^^^ 

5.5182 

-0.0 

1.698 

.9539 

.09312 

.02854 

438.45 

5.4282 

-.188 

1.713 

1.062 

.09324 

.02852 

424.25 

3.3282 

-.568 

1.730 

1.186 

.09337 

.02844 

410.10 

5.2282 

-.520 

1.748 

1.329 

.09352 

.02846 

396.02 

3.1282 

-.648 

1.767 

1.495 

.09369 

.02843 

381.99 

3.0282 

-.755 

1.786 

1.688 

.09389 

.02840 

368.02 

2.9282 

-.844 

1.807 

1.939 

.09411 

.02836 

554.12 

2. 8282/, V 

2.7282'^' 

-.918 

1.828 

2.181 

.09436 

.02832 

340.29 

-.979 

1.851 

2.496 

.09466 

.02827 

326.52 

2.6282 

-1.029 

1.874 

2.872 

.09500 

.02821 

312.85 

2.5282 

-1.069 

1.900 

3.323 

.09539 

.02815 

299.2194 

2.4282 

-1.101 

1.926 

3.868 

.09586 

.02807 

285.69 

2.3282 

-1.127 

1.954 

4.533 

.0^41 

.02799 

272.24 

2.22818 

-1.148 

1.984 

5.350 

.09707 

.02789 

258.89 

2.1282 

-1.164 

2.016 

6.364 

.09787 

.02777 

245.65 

2.0282 

-1.176 

2.050 

7.636 

.09884 

.02764 

252.48 

1.9282 

-1.186 

2.087 

9.252 

.10003 

.02749 

219.43 

1.8282 

-1.193 

2.127 

11.33 

.10152 

.02751 

206.51 

1.7282 

-1.198 

2.171 

14.04 

.10341 

.02710 

193.71 

1.6282 

-1.202 

2.219 

17.65 

.10585 

.02686 

181.05 

1.5282 

-1.205 

2.272 

22.52 

.10904 

.02659 

168.55 

1.42818 

-1.207 

2.332 

29.28 

.11331 

.02628 

156.21 

1.5282 

-1.209 

2.401 

38.91 

.11914 

.02593 

144.08 

1.2282 

-1.210 

2.481 

53.07 

.12728 

.02555 

132.17 

1.1282 

-1.211 

2.577 

74.76 

.13883 

.02513 

120.56 

1.0282 

-1.211 

2.696 

109.5 

.15536 

.02462 

109.58 

.9282 

-1.212 

2.849 

168.3 

.17866 

.02384 

98.85 

.8282 

-1.212 

3.058 

273.2 

.21002 

.02241 

89.36 

.7282 

-1.212 

3.355 

470.8 

.24867 

.01974 

81.30 

.6282 

-1.212 

3.795 

865.8 

.29068 

.01550 

-74.96 

.5282 

-1.212 

4.465 

1717 

.32972 

.01012 

70.40 

.4282 

-1.212 

5.514 

3681 

.35941 

.00478 

67.45 

.5282 

-1.212 

6.968 

7404 

.37363 

.00098 

65. 8r 

.2282 

-1.212 

9.257 

8357 

.36027 

.00001 

65.67 

.1282 

-1. 212 

9.687 

(3)  Free-space  density  P  = 


(1)  See  belov 

(2)  Cavity  radius 
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From  eq.  (4.5),  it  can  be  seen  that  the  {-function  admits  of  am 
arbitrary  mviltiplicative  constant.  This  constant  is  determined  from  the 
known  required  density  of  the  material  behind  the  shock  front  given  by 
the  Ranklne-Hugonlot  conditions.  This  is  given  by  (5.6c).  Prcm  the 
equation  of  state  (6.1), 

P  =  3.5182 

p  ■  384  kb 

Then,  from  (5.10b)  «uid  (5.10c), 

D  -  3.5182  ,  P  «  .02856. 

For  an  initial  cavity  radiiis  of  r  «  I.698  cm,  eq.  (5*8)  gives 
1 

t  ■  r(DP/p)  »  0.868  microseconds. 

This  is  the  value  vhicb  must  be  used  as  the  starting  time  of  the  cavity 
motion  in  order  to  put  the  shock  front  at  the  given  radius.  We  finally 
must  have 

r  -.4 

{  =  »  rt'^  -  (1.698)  (..868  X  10  )  =  451.28. 

11.  Energy  Considerations 

In  the  theory  of  progressing  waves,  an  assunq>tion  of  constant  energy 
was  made  (see  eq.  (4.12)  in  order  to  provide  the  condition  (4.13)  for 
determining  a  and  with  it,  all  the  other  exponents.  The  energy  Integral 
(4.12)  is  extended  between  two  points,  one  of  which  is  located  on  the  shock 
front  {  =  const.  =  end  a  lower  value  1=1^-  The  Integral  path,  such 
as  BC  in  Fig.  3,  may  be  arbitrarily  chosen,  so  long  as  it  terminates  on 


these  two  ctirves. 
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The  energy  in  the  disturbed  part  of  the  solid  will,  however^  change 
with  time  because  its  lower  boundary,  the  ca.vity  surface,  is  not  one  of 
the  family  of  {-curves.  Thus  the  energy  inalues  in  Table  1  represent  an 
integration  of  the  expression  in  (4.12)  taken  along  the  cavity  surface 
curve.  If  we  extend  this  integration  far  enou^  (say  to  100  microsec)  so 
that  point  D  practically  coincides  with  point  C,  the  values  in  the  table 
become  a8yn5)totically  constant,  and  we  have 


From  Table  A  we  see  that  the  energy  does  tend  to  a  constemt  and  we  have 
Just  shown  that  this  limit  is  the  value  of  the  energy  integral  (4,12), 

E  =  1,212  X  10^^  ergs  (11, l) 

If  we  now  suppose  that  all  (or  any  known  fractional  part)  of  the 
energy  given  up  by  the  explosive  is  transmitted  into  the  solid,  then  the 
shock  process  could  be  terminated  when  the  energy  reaches  the  eunount 
available.  It  is  not  possible  to  determine  a  precise  point  of  time  be¬ 
cause  of  the  asymptotic  way  in  which  the  energy  Increases,  It  is  seen^ 
however,  that  E  reaches  90^  of  its  ultimate  value  in  }  microseconds, 
which  is  a  very  short  time  compared  with  the  expansion  process. 

The  post -shock  expansion  presumably  must  take  place  under  constant 
energy  conditions  for  a  "long"  period  of  time,  until  it  is  dissipated  by 
viscosity  of  the  flow,  elastic  waves,  and  other  side  effects. 

For  energy  available  in  the  explosive,  She^r  gives  the  value 
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^otal  “  1“152  k  cal/g 

which,  for  our  explosive  wei^t  of  0.07  lb>  gives 

Oar  calculated  asymptotic  value  (3J..1)  from  the  progressing  wave  integration 
comes  to  9^  of  this.  Thus  we  have  here  an  independent  canparison  to  check 
the  theory. 

12.  Summary  and  Conclxislons 

In  this  paper  we  have  attempted  to  stxidy  the  cavity  fozmation  process 
in  the  metal  hy  determining  how  the  ir^ortant  physical  variables  of  cavity 
radius,  velocity,  pressure,  and  density  vary  with  time  auid  position  near 
the  cavity.  The  most  prominent  geneml  feature  of  the  whole  process  is  the 
short  time  of  the  "shock"  regime  as  compsired  with  the  total  time  of  the  ex¬ 
pansion.  One  general  criterion  of  the  end  of  the  shock  is  when  the  super¬ 
sonic  velocity  of  the  shock  front  drops  to  sonic,  i.e.,  at  the  point  P 
of  Fig.  10,  where  the  slope  attains  the  value  for  elastic  dlsturbemces  in 
the  material.  The  progressing  wave  curve  cannot  be  used  beyon^  this  point 
since  it  would  give  a  subsonic  shock  velocity.  This  situation  has  occurred 
after  0.5  microsec. 

Ve  note  that  the  highest  pressures  and  densities  in  the  metal  are 
located  Just  behind  the  shock  front,  euid  trail  off  with  decreasing  radius 
to  minimum  values  at  the  cavity  boundary.  Ve  note  that  the  equation  of  state 
(6.1)  which  has  been  used  for  the  calculations  has  a  lower  limit  of  p  » 

100  kb.  This  could  also  be  used  as  a  criterion  for  shock  termination  (point 
D,  Pig.  10).  It  is  reached  in  0.7  |x  sec.  These  conditions  thtis  determine 
a  rouj^ily  parallelogram  shaped  region  ODPE  in  the  r,t -plane  for  the 
validity  of  the  progressing  wave  region.  Note  that  the  cavity  has  only 
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ejqianded  O..!.  cm  during  this  period.,  which  is  l/50  of  the  total  ohseirved 
increase  in  radius.-  We  a.re  thus  Justified  in  refe.rri.ng  to  the  shock  process 
as  impulsive,  i..e.,,,  the  later  stages  of  the  prcoess  are  Insensitive  to  many 
features  of  the  shock  partn  Hence  the  progressl.ng  wave  method  of  integra¬ 
tion  remains  valid  for  the  analysis  of  the  shock  zone. 

The  asyn^jtotic  charact-eristics  of  the  progressing  waves  are  thus  not 
of  direct  physical  interest  since  they  do  not  apply  to  the  problem  beyond 
the  region  described  above.  The  expansion  zone,  headed  by  a  wave  travell¬ 
ing  with  the  dilatational  wave  velocity  goes  on  for  at  least  100  p  seconds, 
during  most  of  which  the  metal  continues  to  move  by  fluid  or  plastic  flow. 

The  final  cavity  radius  attained  is  of  great  Interest  to  the  general 
problem  as  this  value  is  directly  observable  on  the  specimens  after  blast. 

In  principle  the  prediction  of  this  radius  should  afford  a  test  of  emy 
theory,  but  the  matter  is  not  so  direct  as  this,  since  several  theories 
are  Involved.  It  is  now  evident  that  the  cavity  formation  process  is  ccm- 
pllcated.  It  starts  under  one  theory  (in  which  the  state  of  the  metal  is 
fSiirly  well  established)  but  terminates  in  a  different  state  of  the  mater¬ 
ial,  ahout  which  information  is  almost  ccmp lately  lacki.ng  Several  mechMi- 
Isois  have  been  suggested  for  terminating  the  cavity  expansion* 

(1)  An  energy -level  criterion 

(2)  A  temperature  criterion 

(3)  A  yield-point  criterion 

Criteria  such  as  (2)  or  (3)  are  tenqptlng  because  they  tend  to  provide 
fairly  definite  marks  as  to  when  the  material  "freezes",  either  when  a  given 
temperature,  or  a  given  pressure  is  reached..  However,  knowledge  of  mater  ¬ 
ials  is  still  too  incomplete  to  solve  this  problem.  The  total  energy  of 
the  moving  material  stops  increasing  after  the  expansion  phase  has  begun, 
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80  there  is  no  change  In  energy.  Furthermore,  any  quantitative  use  of 
energy  halemces  would  require  careful  accounting  of  all  the  energy  loss¬ 
es  as  well.  A  discussion  of  seme  of  these  energy  questions  was  given  in 
the  previous  section. 

Sunmary 

We  summarize  by  noting  that  the  following  four  phenomena  are  coinci¬ 
dent  in  time: 

1.  The  shock-wave  becomes  sonic 

2.  The  pressure  at  the  cavity  surface  drops  to  less  than  100  kb 

3*  The  total  energy  in  the  disturbed  material  reaches  90^  of  Its  maxi¬ 
mum  emd  then  levels  off  aaysptotically. 

4.  The  average  gas  pressure  in  the  cavity  (unlfoni  model)  equals 
that  In  the  metal. 

All  of  these  occur  close  to  2  microseconds  after  initiation  of  the  ex¬ 
plosion.  This  delineates  a  fairly  definite  time  point  of  chsuigeover  of 
conditions.  Xfp  to  this  time  (t  ^  2.5  microsec.  for  the  conditions  of 
this  report)  we  meiy  say  the  effects  of  shock  predominate.  The  progress¬ 
ing  ware  method  furnishes  an  accurate  theory  for  this  regime.  After  this 
time  a  relatively  long  expansion  period  occurs  at  constant  energy  until 
ultimately  temlnated  by  degrsidatlon  processes. 
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APPEHDIX  ~  Numerical  Progr^s 

A.  Calculation  for  the  PjU--dlagrajn 

The  differ entl.al  equation  was  solved  nmerically  on  the  lEM 

707^  at  The  Pennsylvania  State  University  using  a  FORTRAN  program  repro¬ 
duced  below.  A  list  of  the  variables  in  the  notation  of  this  paper,  to¬ 
gether  with  their  symbols  on  the  program  is  given  in  Table  A.  The  method 
used  was  first  to  parametrize  the  curves  P  =  P(u)  in  terms  of  arc  length 
s,  so  as  to  avoid  dl.fficulties  when  the  solution  curves  have  vertical 
tangents-  Thus,  the  element  of  arc  of  the  curve  is  given  by 


ds 


yi  4-  (dP/do)^  dO 


(A-1) 


Substituting  the  expression  ij  8e)  in  this  we  have 
ds 


/ 


!R  +  PS>  =  dU 


',A-2,) 


P^(N  +  PQ}^  *  (R  +  Ps/ 


ds 


yP~(N  +  PQ)^  +  (R  +  PS,'^ 


P<;N  t  PQ’  =  dP 


(A.  5) 


The  square  root  in  the  denominator  Is  never  zero  except  at  the  singular 
points  of  the  differential  equation  To  caic'ilate  the  special  curve  @ 
of  Pig.  7  efficiently  advantage  was  taken  of  the  fact  that  the  source  point 
in  the  rectangle  becomes  a.  "sink"  i  e  a  stable  n-ode  if  the  direction  of 
integration  on  each  curve  is  reversed  Accordingly,  a  starting  point  B 
was  chosen  near  the  saddle  point  ,  -^ith  the  coordinates 

U  =  0  5:990  P  =  C  0.x? 
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(uslng  the  knovm  slope  m  =  -  0.2,  see  Appendix  C,  eq.  (C.2)  and  the 
progmm  I  was  non  with  DS  =  -  0.001.  The  computed  curve  then  ended  very 
accurately  at  the  original  sotirce  point. 

B.  R.T-diagraa  Program 

This  program  is  Independent  of  the  Program  A,  but  actually  incorpor¬ 
ates  the  calculations  of  A  within  itself.  It  calculates  the  "particle 
trajectories"  directly  from  the  initial  conditions.  It  also  calculates 
the  corresponding  P(U)  curve.  For  this  program  it  was  found  best  to 
use  the  density  D  as  Independent  variable.  Incremented  at  unlfom  steps. 
This  is  because  density  is  the  most  sensitive  variable  in  the  portion  of 
the  curves  of  physical  Interest.  (See  Table  1). 


-33- 


TABLE  A  - 

Program  Symbol  Table  for  P, 

U-Dlagram 

Program  Symbol 

Our  Notation  or  Name 

Eq.  Ref. 

Mode 

A 

a 

FLT 

B 

P 

II 

C 

6 

II 

G 

7 

II 

U 

U 

II 

P 

P 

It 

DS 

da 

II 

X 

( 

II 

D 

P 

11 

IMAX 

of  Iterations 

FIX 

EP 

potential  energy 

(8.3) 

FLT 

EK 

kinetic  energy 

(8.5) 

II 

E 

E(t) 

(8.3) 

II 

PN 

N(U) 

(3.8b) 

II 

FQ 

<i(u) 

II 

II 

FR 

R(U) 

It 

II 

FS 

s(u) 

»l 

11 

FT 

P(N+PQ) 

(5.8a) 

II 

IW 

R+PS 

II 

II 

FX 

da/ y  P^(»fPQ)®+(R+PS 

(A.2) 

II 

DP 

dP 

II 

II 

DU 

dU 

(A.3) 

II 

FX 

(U-a)2-7P 

(6.7=) 

It 

DX 

dt 

(6.7a) 

II 

DD 

dp 

(6.5) 

II 

PD 

Pp 

II 

DEP 

(8.3) 

II 

DEK 

II 

ft 
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TABLE  B  -  Additional  Symbols  for  R-T  Diagram 
(See  also  Table  A) 


Program  Symbol 


or  Statement 

Our  Notation  or  Name 

Mode 

DO 

*>o 

FLT 

R 

R 

tt 

St.  47 

Initial  U 

(11.5a) 

tl 

St.  48 

Initial  P 

(11.3c) 

It 

St.  54 

Initial  p 

(11.5b) 

II 

EX 

1/a 

It 

T 

t 

II 

V 

u 

II 

PRES 

P 

II 

DT 

dt 

II 
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C.  Singular  Points  of  the  P.U-Plane 
The  differential  equation  (3* 8b.) 

dP  P[M(U)  +  pa(u)l 
dU  '  r(U)  +  l4(U) 

is  said  to  have  a  singular  point  whenever  ^  =  g'  ^  Using  the  expressions 
(3- 8b),  this  condition  occurs  If 

(51)  P  a  0,  U  =  0 

(52)  P  a  0,  U  a  a 

(53)  P  a  0,  U  a  1 

and  where 


This  condition  leads  to  the  cubic  equation 

+  AgU^  +  Aj^U  +  A^  -  0 

where 

Aj  -  27  (1  -  37) 

Ag  -  37^(3a  -  1)  +  7  (7  -  5a)  -  r(3B  +  4p)  +  B  +  2p 
Aj^  -  7(308  +  6c3p  -  2p  -  lo)  +  (2  -  a)  (B  +  2p) 

A^  -  -  aa  (b  +  2P) 

For  the  partlcxilar  nunerlcal  constants  given  by  (5.9)  suid  7  »  7>6,  the 
point  S4  is  given  by 

P  a  0.02884  ,  U  a  0.0916 

For  the  given  constsuits  and  7  a  7.6  for  alxsnlnuin,  the  renalnlng  roots 
of  the  cubic  are  imaginary. 
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Near  any  isolated  singular  point  P^)  we  nay  approximate  the 

differential  equation  by  a  linear  fractional  expression  of  the  form 


dP  ^  a  X  +  by 
dU  c  X  +  dy 

where  x  »  U  -  y  ■  P  "  Pq*  See  Stoker  [13]>  P*  36-45,  where  tests 
are  given  for  determining  the  type  of  singular  point.  A  summary  of  the 
resiilts  is  as  follows: 


Singular  Point 

U  ■  0,  P  «  0 
0  ’  o 


U  -  a,  P  -  0 
o  ’  o 

U  -  1.0,  P  -  0 
0  '  o 


Coefficients 

a  =  0,  b  »  -  20 
c  -  -  O,  O  -  2P 


Classification 
Stable  Node 


(see  below)  Saddle  Point 

a  »  0,  +  2r  Stable  Node 

c  =  0-1,  d  «  37  +  2P 


U  -  .0916,  P^  .  .02884 
0  o 


Uhstable  Node 


Since  the  point  «  O,  ^  special  importance  (end-point  of 

the  solution  curve)  it  was  examined  directly.  Near  U  =  O  we  have 
the  power-series  ejqpeuision. 

Dividing  both  numerator  suid  denominator  in  the  differential  equation  (4.8) 
by  U  -  o  and  putting 


m  «  P/(U  -  o) 


(C.2) 


we  nay  write 
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^  _  FnCu)  +  m  (26  +  2y(U  -  a)1 

dU  “  U(i  -  U)  +  (2P  +  5U7)in 


(C.5) 


We  new  inquire  whether  a  solution  curve  can  approach  the  point  U  *  a 
with  a  definite  limiting  slope  (i.e.  be  locally  a  straight-line).  Letting 
U  -^a  in  (C.5)  putting  this  expression  in  (C.l),  emd  dividing  by  (U  -  a) 
we  obtain  the  equation  for  m: 


m 


m 


a(a  -  l)(y  -  1)  +  2g  m 
a  (1  -a]  +  m(2p  +  3ra) 


This  has  the  following  three  solutions: 

1)  m  =  09 

2)  m  s  0 

3)  m  »  (a  -  l)/3  for  ^  =  a  -  1 

Thusthta*  are  three  solution  CTirves  which  pass  through  the  point.  The 
first  two  are  the  straight  lines  U  »  a  and  P  >  0  respectively,  which 
are  singular  solutions  to  the  differential  equation.  The  roaalning  slope 
is  m  s  -  0.2  for  a  s  0.4^.  This  is  the  limiting  slope  of  the  solution 
curve  0  of  Fig.  7-  In  the  same  figure,  all  the  solution  curves  between 
(0  and  the  horizontal  line  P  =  0  end  up  at  the  origin,  while  those  be¬ 
tween  (0  and  the  vertical  line  U  »  a,  go  to  infinity  aisynptotlc  to  the 


latter . 
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COMPILE  RUN  FORTRAN 
BRL-PENN  STATE  PROJ  BLAST  WAVES 
P-U  DIAGRAM  PROGRAM 
DAVIDS-CALVIT 
DIMENSION  IOENT(8) 

1  READ  600*A*B.C»G* IDENT 
600  FORMAT(4F10.0*8AS  I 

IF(AI2 •30<2 

2  READ  801«UtPt0S«X*D«IMAX 
801  FORMAT(5F10.0.I10) 

200  READ  807.EP.EK.E 
807  FORMAT ( 3F10.0 » 

PRINT  605*I0ENT 

805  FORMAT (1H1,5X.31HBRL-PENN  STATE  PROJ  BLAST  WAVES»4Xt8A5) 
PRINT  806 

806  FORMAT! IHOtlH I «6XtlHU*9X»lHP.9X»lHX*9X*lHD*14X*2HEP« 
113X«2HEX*13X*1HE«9X*2HP0) 

3  DO  26I>1«IMAX 

4  FN-6*U*(3.0»A-1«0-2.0»U  »♦ ( 3. 0-A > *U-2 .0*A 

5  FQ»(2.0*B-(G-1.0)«C»/(U-A)+2.0*G 

6  FR-U*(U-A)*( 1 .0-U» 

7  FS«C+2.0*B+3.0*U«G 

8  FT«P*( FN+P*FOI 

9  FW«FR+P*FS 

10  FX«DS/SQRTF(FT*FT+FW*FWI 

11  OP-FT*FX 

12  OU-FW*FX 

13  FY»(U-A>*(U-A|-G*P 

14  0X-X*FY*FX 

15  00«-0*(FW+3.0*U*FY>*FX/(U-A» 

P0»P*0 

16  PRINT  803tl .U.P»X.OfEP,EKtE»PD 

803  FORMATdH  ♦  I  3  ♦4F  10.5  »  3F 1 5. 7 , 1 F 10. 5  ) 

17  U-U+OU 

18  IF(U)27,27,19 

19  P-P+OP 

20  X-X«OX 

21  D>0+00 

22  0EP«4.0*3.141593*P*0*!X**4)#0X/(G-1*0) 

23  OEK»4.0*3.141593* (U*U)*O*!X*»4)*0X/2.0 

24  EP-EP+OEP 

25  EK-EK+OEK 

26  E-EP+EK 

27  GO  TO  1 

30  STOP 

31  END 


n  r»  r» 
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COMPILE  RUN  FORTRAN 
BRL-PENN  STATE  PROJ  BLAST  WAVES 
R-T  CALCULATION  PROGRAM 
DAVIOS-CALVIT 

53  DIMENSION  I  DENT ( 6 ) • JOENT ( 6 ) 

1  READ  a00*A»B.C.GtDOi lOENT 
800  FORMAT(5F10.0t6A5 ) 

2  READ  809»X,R,D0» IMAX. JOENT 

809  FORMAT (lF10.0t2F10.0. I 10»8A5) 

43  IF( IMAX)30»30.44 

44  PRINT  805 

805  FORMAT IlHl ♦5X»31HBRL-PENN  STATE  PROJ  BLAST  WAVESI 

45  PRINT  810iA.B»C*GfDO*I0ENT 

810  FORMAT  I lH0t5X,5F15.5.5X.6A5) 

46  PRINT  813»X«R«DOtIMAXtJOENT 

813  FORMAT! lH0t5X«3F15.5iI10*5X»8A5) 

47  U«A*( 1*0-(G-1.0I/(G+1.0) I 

48  P- IA*A*IG-1.0)/ (G+1.0) ) « ( 1 .0- (G-1 .0 )/ (6+1.0 ) » 

54  D-(6+1.0)/(G-1.0J*DO 

49  EP-0 

50  EX-0 

51  E-0 

52  PRINT  806 

806  FORMAT!  IHOtlX.lHI tTXtlHT  tllX*lHR»12X»lHV»llX*7HPRES  XB*8X » IHD . 1 1X« 
llHU«llX«lHPillX»lHX«12XtlHE| 

3  00  260  I-ltIMAX 

4  FN<»G*U*!3.0*A-1.0-2.0*U)+!3.0-AI*U-2.0*A 

5  FQ-|2.0*B-!G-1.0»*Cl/!U-A)+2.0*6 

6  FR«U*(U-A>*! 1,0-U> 

7  FS-C+2.0*B+3.0*U*G 

8  FT-P*!FN+P*FQl 

9  FWFR+P+FS 

13  FY«!U-A)*!U-A|-G»P 

900  FX--00*!U-A»/!0*!FW+3.0*U*FY) ) 

11  DP»FT*FX 

12  OU*FW*FX 

33  EX-1. 0/A 
330  T»!R/XI**EX 

14  DX»X*FY*FX 

34  V-R+U/T 

37  PRES«!R/T 1**2 *P*0»! 10.0** 1-9 U 
160  PRINT  808 . 1  tT.R.V.PRESfO.U.PtX.E 

808  FORMATIIH  . I  3 » 1 E 12 . 4 . IF  12 • 7. lE 14 .5 . 1 F 14. 5 .4F 12 • 5 . 1 E15 . 5 ) 

340  DT-DX/ ! (-A/T+V/R) *X ) 

19  IF!PI27.35,35 

35  T-T+DT 

36  R-R+V*DT 

17  U-U+OU 

18  P-P+DP 

20  X-X+DX 

21  D-O+DO 

22  OEP=4.0*3.141593*P*D*!X**4|*OX/(G-1 .0) 

23  DEK-4.0*3.141593* !U*U>*0*!X«»4)*DX/2.0 

24  EP-EP+OEP 

25  EX  =  EX  +  OEIC 

26  E-EP+EX 
260  CONTINUE 

27  GO  TO  2 

30  STOP 

31  END 


FIG.  I  “  ALUMINUM  SPHERE  -  PRE 


reflected  wave 
Zone  B  -  Heavily  damoged 
Zone  C  —  "Fluid”  zone 
Note  -  Oimenslont  ore  approximate 


FIG.  2  -  ALUMINUM  SPHERE  -  POST  SHOT. 
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FIG  3  -  CAVITY  EXPANSION  AND  PARTICLE  TRAJECTORIES 
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4  -  EQUATION  OF  STATE  OF 
ALUMINUM 


400 
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FIG  5 -EQUATION  OF  STATE,  EXPLOSION  PRODUCTS  OF 
PENTOLITE  (R.  E.  SHEAR) 


Gas-Aluminum  Inlerface 
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/>•/>. 


R  (t) 

(a)  Uniform  Density 


(b)  Nonuniform.  Detonation  Wove  Remains  in 
Contact  With  Metal  Boundary 


(c)  Detonation  Wove  Reflected  froth  Metal 
Boundary 


FIG.  6  -  POSSIBLE  MODELS  FOR  DENSITY  VARIATION 
OF  GAS  IN  EXPANDING  CAVITY 
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-P,U  DIAGRAM.  PROGRESSING  WAVES  :  AL .  X  =  7.6 


BEHAVIOR  OF  SOLUTION  CURVES  NEAR  U=  0916,  P  =.02884  foe 
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fig.  9 


(a)  FINITE  SHOCK  STRENGTH 


(b)  INFINITE  SHOCK  STRENGTH 


-  POSTULATED  PARTICLE  PATHS 
ACROSS  SHOCK  FRONTS. 
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DIAGRAM  OF  SHOCK  REGION  FOR  ALUMINUM  SPHERE 


R-T  DIAGRAM  FOR  CAVITY  SURFACE  FOR  DIFFERENT 
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(q>l) 
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FIG.  12  -  PRESSURE  DECAY  ON  CAVITY  SURFACE  IN  AL 
y  •  7.6  ,  P  -  2.70  gm/cm* 
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FIG  14  -  VARIATION  OF  DENSITY  AT  CAVITY  SURFACE  ; 
AL  :  y  =  7.6  Ao  =  2.70  gm/cm* 


